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ABSTRACT 

The rapidly oscillating Ap (roAp) star 10 Aql shows one of the lowest photometric pulsation 
amplitudes and is characterized by an unusual spectroscopic pulsational behavior compared to 
other roAp stars. In summer 2006 this star became target of an intense observing campaign, 
that combined ground-based spectroscopy with space photometry obtained with the MOST 
satellite. More than 1000 spectra were taken during 7 nights over a time span of 21 days 
with high-resolution spectrographs at the 8-m ESO VLT and 3.6-m TNG telescopes giving 
access to radial velocity variations of about 150 lines from different chemical species. A 
comparison of pulsation signatures in lines formed at different atmospheric heights allowed 
us to resolve the vertical structure of individual pulsation modes in 10 Aql which is the first 
time for a multiperiodic roAp star. Taking advantage of the clear oscillation patterns seen in 
a number of rare-earth ions and using the contemporaneous MOST photometry to resolve 
aliasing in the radial velocity measurements, we improve also the determination of pulsation 
frequencies. The inferred propagation of pulsation waves in 10 Aql is qualitatively similar 
to other roAp stars: pulsation amplitudes become measurable in the layers where Y and Eu 
are concentrated, increase in layers where the Ha core is formed, reach a maximum of 200- 
300m s -1 in the layers probed by Ce, Sm, Dy lines and then decrease to 20-50ms~ 1 in 
the layers where Ndlll and Prill lines are formed. A unique pulsation feature of 10 Aql is a 
second pulsation maximum indicated by Tb III lines which form in the uppermost atmospheric 
layers and oscillate with amplitudes of up to 350 ms -1 . The dramatic decline of pulsations in 
the atmospheric layers probed by the strong Pr III and Nd III lines accounts for the apparent 
peculiarity of 10 Aql when compared to other roAp stars. The phase-amplitude diagrams and 
bisector measurements of the Nd III 5 102 A line reveal a rapid change of phase and amplitude 
with height for all three main pulsation modes, indicating the presence of a pulsation node in 
the stellar atmosphere. Finally, we report the discovery of a puzzling asymmetry of the strong 
Ndlll lines with their blue wing extending up to — 50kms _1 , which is about 25 times the 
estimated value of v e sin i. 

Key words: stars: atmospheres - stars: chemically peculiar - stars: magnetic fields - stars: 
oscillations - stars: individual: 10 Aql 



1 INTRODUCTION 

More than 30 members among the group of late A magnetic 
chemically peculiar stars exhibit high-overtone, low-degree, non- 
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radial p-mode pulsations w ith periods in the range of 6-21 min- 
utes I IKurtz & Martinez 2000). These rapidly oscillating Ap (roAp) 
stars are characterized by strong global magnetic fields with a po- 
lar strengths of the order of 1-10 kG. The atmospheres of roAp 
stars are enriched with heavy elements, brought up from the stel- 
lar interior by diffusion. Conspicuous lateral and vertical variations 
of chemical abundances in Ap stars and the prominence of the 
magnetic field signatures make them primary targets for detailed 
investigations of the magnetic topology and magnetically driven 
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formation of struc tures in stellar atmospheres dKochukhovll2004al ; 
lRvabchikovj|2004h . 

The observed pulsational amplitudes of roAp stars are mod- 
ulated according to the visible magnetic fiel d structure. T his ob- 
servation led to the oblique pulsator model dKurtz|[l982h . where 
axisymmetric 1 = 1 modes are aligned with the magnetic field 
axis, whi ch itself is oblique to the axis of stellar rotation. Calcu- 
lations bv lSaio & Gautschvl f2004h showed that dipolar modes are 
significantly distorted by a magnetic field of kG-strength, but retain 
their axisymmetric character. On the other hand, the structure of 
pulsational perturbations can inc lude significant non-axisymm etric 
components in weak-field stars (Bigot & Dziembowski 2002]). The 
indirect pulsation Dop pler imaging of osc illations in the prototype 
roAp stars HR3831 jKochukhovlboO^ vindicated the oblique 
pulsator model and for the first time characterized observationally 
the magnetic distortion of the global p-modes. 

The roAp stars are key objects for asteroseismology, which 
presently is the most powerful tool for testing theories of stellar 
structure and evolution. The classical asteroseismic analysis, uti- 
lizing precise frequency measurements, helps to constrain the lu- 
minosity and internal chemical composition of pulsating magnetic 
stars (e.g. Matth ews, Kurtz & Martinez l 19991 ; Cunha, Fernandes & 
Monteiro 120031) . Recent time-resolved spectroscopic observations 
of roA p stars uitilizing large telescopes (see review by Kochukhov 
l2007al) demonstrated the possibility of another type of asteroseis- 
mic investigation, which is focused on the upper atmospheric lay- 
ers. A significant chemical stratification and a short vertical wave- 
length of pulsation modes lead to a remarkable diversity of the pul- 
sational characteristics observed in individual spectral lines, no- 
tably of absorption lines of the rare-earth elements (REE) (e.g., 



Kochukhov & Rvabchikova 2001a; Mkrtichian, Hatzes & Kanaan 
2003; Rvabchikova et al. 2007a). The information extracted from 



the lines formed at different optical depths opens access to different 
modes and can be combined to yi eld a vertical tomographic map of 
the pulsating stellar atmosphere (Rvabchikova et al. 2007b). Thus, 
applying Doppler imaging techniques and pulsation tomography to 
line profile variation in roAp stars pro mises an unpreced ented 3-D 
picture of magnetoacoustic pulsations dKochukhovll2005t) . 

These spectacular results have been obtained despite the fact 
that most spectroscopic observations of roAp stars were obtained 
in a snapshot mode and hence cover only a few hours of stel- 
lar oscillations. A detailed frequency analysis is highly ambiguous 
with such limited data, because of the well known aliasing prob- 
lem. It can only marginally improved by combining observations 
from several nights, but which rarely are available when large tele- 
scopes had to be used. Therefore, despite an outstanding clarity of 
the pulsation curves often produced with spectroscopy, for most 
multiperiodic roAp stars they correspond to an unresolved mix- 
ture of different modes and thus are exceedingly difficult to inter- 
pret. This fundamental problem can be alleviated by simultaneous 
photometric observations with a high duty cycle, as it was suc- 
cessfully demonstrated by a combination of space photometry of 
the roAp star HP 247 12 with contem poraneous VLT/UVES spec- 
troscopy (Rvabchikova et al. 2007a). In this paper we present an 
analysis of a much more extensive spectroscopic time series for 
another pulsating Ap star, 10 Aql, which was observed simultane- 
ously from ground and space. 

lOAql (HR7167 HP 176232, HIP 9 3179) is one of the 
brightest roAp stars. iRvabchikova et all (2000) performed its 
model atmosphere and abundance analysis using high-resolution 
spectra. They have estimated T c g = 7550 K, log g = 4.0, and de- 
rived M = 2.0 ± 0.2M Q , R = 2.5 ± 0.2 R Q from the compari- 



son of the stellar parameters with theoretical evolutionary tracks. 
iKochukhov & Bagnulol d2006t) have included 10 Aql in their study 
of the evolutionary state of magnetic chemically peculiar stars. 
Adopting a photometric temperature T c ff = 7900 K, they found 
Af=1.95 ± 0.04 M Q , log L = 1.32 ± 0.05 and a fractional age 
of 64-76% of the main sequence lifetime. The abundance pattern 
of 10 Aql is characteristic of a cool Ap star. A notable spectral 
anomaly is the overabundance of the doubly ionized REEs, Prill 
and Ndm. 

A long itudinal magne t ic field of abou t 50 G was detected in 
10 Aq l bv lBabcockl jl958h . |Prestonl jl970h and IRvabchikova et all 
(2005). Using m agnetic Zeeman splitting in unpolarized spectra 



IKochukhov et al.l d2002l) provided the first measurement of the 
mean magnetic field modulus in 10 Aql, (73) = 1 .5 ± .1 kG, which 
later was confirmed by Leone, Vacca & Stift J2003I) . and he es- 
timated the projected rotational velocity v c sin i = 2.0 km s" 1 . The 
sharpness of spectral lines and the absence of measurable variations 
of the mean longitudinal magnetic field indicate that 10 Aql has a 

long rotational period. 

Recently, Ryabchikova, Kochukhov & Bagnulo d2008l) have 
studied the vertical stratification and isotope anomaly of Ca in 
a number of Ap stars, including 10 Aql. These authors found a 
concentration of Ca in deeper atmospheric layers and invoked a 
vertical separation of the Ca isotopes to explain the shape of the 
Call infrared triplet lines. A comprehensive chemical stratifica- 
tion analysis of 10 Aql is currently underway by Nesvacil et al (in 
preparation), but already a preliminary result (Nesvacil, Weiss & 
Kochukhov 2008) shows that Si, Ca, Cr, Fe, and Sr share a quali- 
tatively similar vertical abundance distribution, characterized by a 
rapid increase of the element abundances below logrsooo ~ — 1. 
A re-determination of the atmospheric parameters taking into ac- 
count individual chemical abundances and str atification supports 
the sp ectroscopic T c g ~ 7600 K obtained by IRvabchikova et all 
J2000h - 

Photometric variations with a period close to 11.4min and 
an amplitude be l ow . 5mmag were dis c overe d in 10 Aql by 
iHeller & Kramed dl988t> . IHeller & Krameil dl990h identified three 
pulsation modes with periods of 1 1 .06-13.45 min, based on 26 h of 
high-speed p hotometry in the /3-band . Pespite an aliasing problem 
in their data, IHeller & Kramed d 1990S) were able to determine the 
large frequency separa tion, A y = 50.6 fjHz. Belmonte, Martinez 
Roger & Roca Cortes dl99ll) collected 47 h of time-resolved pho- 
tometry in the J-band, confirming oscillations in 10 Aql. Moreover, 
their study suggested that the amplitude of the infrared photomet- 
ric variations significantly exceeds variations in the 73-band, which 
is unusual for a roAp star (e .g. Martinez, Sekiguchi & Hashimoto 
1 19941; iMatthews et ai]|l996t) . However, no oscillations were de- 
tected in the follow-up photometric monitoring conducted in the 
//-band (Belmonte, Kreidl & Martinez Roger ll992j) . 

Radial velocity variations with pre viously known os c illatio n 
periods were discovered in 10 Aql by Kochu khov et al.l d2002h . 
Over two nights they observed the star during a total of 8 h in a 
short wavelength region centred at A 6150 A. Unlike the majority 
of other roAp stars, 10 Aql showed only a weak variability with 
an amplitude of ~30ms _1 in the strong Ndm 6145 A line, but 
exhibited amplitudes at the level of 80-130 ms _1 in the lines of 
first REE ions, Gd II and Eu II. Aiming to extend t he stu dy of pul- 
sational variation of 10 Aql. IHatzes & Mkrtichian! d2005l) analysed 
spectra covering 1 1.6 h of time-resolved observations during three 
nights. Somewh at surprisingly, despite a bro ad wavelength cover- 
age of their data, Hatzes & Mkrtichian (2005) were able to confirm 
the presence of RV variations in only 5 individual spectral lines, 
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Table 1. Pulsa tion frequencies o f 10 Aql identified in the MOST photomet- 
ric campaign iHuberetal. 2008). The numbers in brackets give the esti- 
mated error in units of the last significant digit. The two candidate frequen- 
cies, fi and fj, are possibly present in the MOST data but could not be 
definitely confirmed. 



id 


v (mHz) 


P (min) 


A (mmag) 


h 


1.44786(3) 


11.5112(2) 


0.17(1) 


h 


1.39691(3) 


11.9311(3) 


0.15(1) 


h 


1.42709(4) 


11.6788(3) 


0.12(1) 


h 


1.3662(1) 


12.199(1) 


0.03(2) 


h 


1.4686(1) 


11.349(1) 


0.03(2) 



two of which remained unidentified. The authors found the highest 
RV amplitude of 400 ms" 1 in an unidentified line at A 5471.42 A 
but could not detect pulsations in strong doubly ionized REE lines 
which show prominent variability in other roAp stars. In addition 
to its remarkably low-amplitude photometric pulsational variations, 
10 Aql seems to exhibit an unusual spectroscopic pulsational be- 
haviour, making this star perhaps an intermediate object between 
high-amplitude roAp stars and apparently constant non-pulsating 
magnetic Ap stars. 

10 Aql was chosen as a target for an extensive observing 
campaign by the C anadian photometric space telescope MOST 
(Walker et al. 2003). The star was observed as a primary target in 
June- July 2006, and nearly 120000 10 s exposures were collected 
du ring 3 1 d at a samp ling interval of 20 s. The frequency analysis 
bv lHuber et al.l |2008) revealed three definite periodicities and two 
candidate frequencies (see Table[T](. Hube r et al. showed that the 
two largest amplitude peaks reported by iHeller & Kramer! Jl990h 
are Id -1 aliases of the intrinsic stellar frequencies /i and /2. The 
third significant MOST frequency, f$ — 1.43 mHz, has not been 
detected in previous photometric observations of 10 Aql. The lack 
of rotational modulation signal in the MOST data confirms the sus- 
pected slow rotation of the star. iHuber et al. (2008) derived an im- 
proved value of the large frequency separation, Av — 50.95 /^Hz, 
and have compared the observed frequencies with the predictions 
of the theo retical mo dels of non-adiabatic, magnetically distorted 
p-modes dSaid2005h . 

With the aim to improve the frequency analysis, study line 
profile variations and perform a tomographic analysis of the pul- 
sations in the atmosphere of 10 Aql, we have organized a ground- 
based spectroscopic observing campaign simultaneously with the 
MOST observations. Using high-resolution spectrographs at 4- and 
8-m telescopes, we have collected a superb observational material 
with very high signal-to-noise ratio, wide wavelength coverage and 
a high spectral resolving power. This combination of the data qual- 
ity and time coverage has been never achieved before for any roAp 
star. 

Preli minary results of our pulsation study of 10 Aql were re- 
ported bv lSachkov et al.l ( 120071) . In that paper we outlined the fre- 
quency analysis procedure and discussed the main pulsation prop- 
erties of 10 Aql, such as the presence of a node-like behaviour in the 
bisector variations, phases and amplitudes of radial velocity varia- 
tion fo r different REE ions. Subsequently, Elkin, Kurtz & Mathys 
(2008) published an analysis of a short spectroscopic UVES time 
series of 10 Aql. These data, covering only 2 hours, were analysed 
without taking the multiperiodic character into account. The pulsa- 
tion frequency of 1.428 mHz adopted by Elkin et al. coincides with 
neither of the two highest-amplitude pulsation modes in our data, 



but is close to the frequency fa identifyed by |HuberetaIU2008h 
and by our spectroscopic analysis (see below). 

The present paper is organized as follows. In Sect.|2]we de- 
scribe spectroscopic observations and data reduction. Section[3]de- 
scribes methodology of the line identification and radial velocity 
measurements. The frequency analysis is outlined in Sect.|4] We es- 
tablish a relation between photometric and spectroscopic variations 
of 10 Aql in Sect.[5] Line profile variations, indicating a pulsation 
node in the stellar atmosphere, and the discovery of a remarkable 
asymmetry in the Nd III lines are presented in SectJJj] The paper 
ends with conclusions and discussion in Sect.[9] 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Time-series of UVES spectra 

The main set of observational material analysed in this paper was 
obtained at the ESO 8-m VLT UT2 telescope using the UVES spec- 
trograph. 10 Aql was observed for about 4 h during each of the fol- 
lowing four nights in 2006: July 3, 9, 15, and 17. Each subset of ob- 
servations consisted of 211 stellar exposures, followed by a ThAr 
lamp spectrum. We have used UVES in the 600 nm red-arm set- 
ting, that provided access to the 4980-6980 A spectral window. The 
wavelength coverage is complete, except for a ~ 100 A gap cen- 
tred at A 6000 A. To facilitate the analysis of the narrow line pro- 
files of 10 Aql and to improve the wavelength stability we used the 
UVES image slicer No. 3, resulting in a spectral resolving power 
of A/ A A « 115 000. 

10 Aql was observed with exposure times of 40 s. The fast 
readout mode of the UVES CCDs (625 kpix s — 1 , 4-port, low gain) 
reduces the overhead to about 29 s, resulting in a time resolution 
of 69 s. The typical signal-to-noise ratio of individual spectra is 
250-300 at A = 5000-5500 A. 

Reduction of the UVES spectra utilized a set of echelle spec- 
tral pr ocessing routines developed by Tsymbal, Lyashko, & Weiss 
(2003) and Lyashko, Tsymbal & Makaganuik 120071) . Their code 
performs the standard reduction steps (averaging of calibration 
frames, determination of the echelle order position, bias subtrac- 
tion, flat fielding, extraction of 1-D spectra and wavelength cal- 
ibration) and is optimized to treat the UVES slicer spectra. Fur- 
ther details of the application of this spectral reduction package 
to the time-resolved UVES observations of roAp stars are given 
bv lRvabchikova e t al. (2007b). This paper also details the final re- 
duction step in which we improve continuum normalization of the 
individual time-resolved spectra obtained in the same night. Obser- 
vations were obtained in different nights, hence we have carried out 
an additional continuum rectification correction to achieve a global 
consistency of the continua in all 844 UVES spectra. 

Our dataset of the UVES time-resolved spectra was comple- 
mented with the observations obtained in the context of program 
077.D-0150 dElkin et alj|2008l) on July 24, which were taken with 
the instrumental setup identical to ours and consist of 105 expo- 
sures obtained with the same exposure time and sampling as in 
our observations. These spectra along with the appropriate calibra- 
tions were retrieved from the ESO Archive and were subsequently 
processed in the same way as our own observations. Thus, a total 
of 949 high-quality UVES spectra are available for the pulsational 
analysis of 10 Aql. 
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Table 2. Journal of spectroscopic observations of 10 Aql. 



Date 


Start HJD 


End HJD 


Spectral range 


N spec. 


Exposure 


Time 


Typical 


Instr. 




(245 0000+) 


(245 0000+) 


(A) 




time (s) 


resolution (s) 


S/N 




03/07/2006 


3919.60523 


3919.77510 


4980-6980 


211 


40 


69 


330 


UVES 


09/07/2006 


3925.64884 


3925.81471 


4980-6980 


211 


40 


69 


260 


UVES 


15/07/2006 


3931.62381 


3931.79193 


4980-6980 


211 


40 


69 


240 


UVES 


17/07/2006 


3933.63894 


3933.80793 


4980-6980 


211 


40 


69 


300 


UVES 


24/07/2006 


3940.64279 


3940.72682 


4980-6980 


105 


40 


69 


300 


UVES 


14/07/2006 


3930.52693 


3930.61096 


4572-7922 


59 


70 


120 


150 


SARG 


15/07/2006 


3931.52970 


3931.66690 


4572-7922 


86 


70 


120 


90 


SARG 


16/07/2006 


3932.49220 


3932.58039 


4572-7922 


61 


70 


120 


140 


SARG 



2.2 Time-series of SARG spectra 

In addition to the UVES time series, 10 Aql was observed at high 
resolution with the spectrograph SARG at the 3.55-m Telescopio 
Nazionale Galileo (TNG) at the Observatorio del Roque de los 
Muchachos (La Palma, Spain) in the three consecutive nights of 
July 14-16, 2006. The total number of echelle spectra obtained 
with the SARG is 206. The star was monitored for 2.0-3.3 h using 
70 s exposure times in each of the three nights. Including the 50 s 
overhead, the time resolution of these observations is 120 s. The 
SARG spectra cover the range of 4572-7922 A without gaps with 
a resolving power of about 86 000. The typical S/N ratio achieved 
for individual exposures is 90-150 at A = 5000-5500 A. 

The spectra were reduced with standard procedures for spec- 
troscopic observations which are part of the Image Reduction and 
Analysis Fac ility package of NAOA in t he same way as in Leone 
& Catanzaro iLeone & Catanzarcl ( 120041) . The post-processing of 
1-D extracted SARG spectra was performed consistently with the 
UVES spectra. 



3 LINE IDENTIFICATION AND RADIAL VELOCITY 
MEASUREMENTS 

To perform a careful line identification and to choose suitable 
lines for pulsation analyses, we have synthesised the whole spec- 
tral region with the model atmos phere parameters T e g = 7 550 K, 
log g = 4.0, and abundances from Rvabchi kova et al. i 2000h. The 
mean magnetic field modulus of (B) = 1.5 kG iKochukhov et al.l 
l2002h wa s adopted. Atomic data were extracted from the VALD 



database dKupka et alj | 1999|). the revised atom ic parameters for 
Nd III were taken from Rvabchik ova et al] {2006), and the synthetic 
spectrum calculations were carried out with the S YNTHMAG code 
( |Piskunov|[l999l ; lKochukhovll2007bh . 

We have measured ra dial velocities using the centr e -of- gravity 
technique as outlined by iKochukhov & Rvabchikov3 d200 1 ah . A 
relatively low pulsation amplitude and a small number of spectral 
lines showing pulsation signatures seriously limited the usability 
of identification lists developed for other roAp stars. Therefore, we 
started by measuring practically all spectral lines (about 2000) in 
the July 3 time series in order to find lines with useful pulsation 
signatures. It turns out that the total number of such lines, both 
blended and unblended, is about 150. Most of the other lines belong 
to iron-peak elements which do not show significant RV variations 
or are blended to an unacceptable degree. 

It is a well-known fact that spectroscopic pulsational vari- 
ability of typical roAp stars is dominated by the lines of rare- 
earth ions, especially those of singly and doubly ionized Pr 



and Nd, which are strong and numerous in the roAp s pectra 
(see, for example, Savanov, Malan ushenko & Rvabchikova 1 1 9991 : 
IKochu khov & Rvabc hikovall200l3) . 10 Aql is different in this re- 
spect from other roAp stars. REE lines are weak, in particular those 
of the first ions. Their equivalent widths are confined between 2 
and 12 mA. The Pr and Nd overabundance inferred from the lines 
of the first ions does not exceed 1.0 dex. Radial velocity measure- 
ments from these very weak lines often have insufficient precision 
for the purpose of our investigation. Therefore, we used for our fur- 
ther analyses only the lines with an equivalent width larger than 
5 mA. The line identification list compiled for the July 3 spec- 
tra was then applied to all other UVES observations, using con- 
sistent central wavelengths and measurement windows. For strong 
and intermediate-strong lines our measurements of individual ra- 
dial velocities have an internal accuracy of up to 10 ms _1 . 

We have unambiguously detected pulsations in the core of Ha 
and in the lines of Y II, La II, Ce II, Pr III, Nd II, Nd III, Sm II, Eu II, 
Gdll, Tbm, and Dy III. All lines show strong changes in RV am- 
plitude, indicative of beating between several excited modes. Typ- 
ical observed semi-amplitudes, averaged over the whole duration 
of our spectroscopic monitoring, range from 50 to 150 m s _ . Dur- 
ing the periods of constructive interference amplitudes approach 
400 m s _1 for some REE lines. 

Remarkably large RV amplitudes - 350 ms' 1 on average - 
were measured for the Dy III 5345, 5556, 5730 A lines and for the 
two unidentified lines at A 5373.02 and 5471.40 A. The identifi- 
cation of 5345 and 5556 A lines is based on the unpublished data 
kindly provided by A. Ryabtsev (Institute of Spectroscopy RAS). 
All these lines have equivalent widths in the range of 9-12 mA 
and their depths are less than 9% of the continuum. The pulsational 
characteristics of the two unidentified features are similar to those 
of Dy III and a few Ce II lines. Without any doubt both unidenti- 
fied lines be long to the REE species. B oth features are measured in 
HD 24712 (Rva bchikova et al.| [2007a) and their pulsational char- 
acteristics are much closer to DyHI lines than to Cell lines. In 
7 Equ, too, Dy III 5730 A and 5373.02, 5471.40 A lines have sim- 
ilar pulsational amplitudes and phases, which differ significantly 
from those of Cell lines. Thus, we suggest that the unidentified 
lines at 5373.02 and 5471.40 A5471.40 belong to the Dy III spec- 
trum. 

The strong Ndlll and Prill lines in 10 Aql, that usually ex- 
hibit the largest amplitudes in roAp stars, have amplitudes below 
100 ms _1 and probably form in the nodal area (see Sect. [6}. The 
strongest Ndll lines in the observed spectral range have equiva- 
lent widths smaller than 7 mA, while PrII lines are not measur- 
able at all. Hence, the blending problem in this star becomes more 
serious. Even slight blending of a REE feature may sign ificantly 
decrease the inferred amplitude iKochukhov et al.l 2002). Taking 
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into account the observed pulsation amplitudes and accuracy of RV 
measurements (defined primarily by the spectral line strength and 
central wavelength), we have selected 65 lines for the frequency 
analysis and the study of the pulsational wave propagation through 
the stellar atmosphere. This line list is provided in Table|4] 

We also report on the analysis of a subset of non-pulsating 
lines belonging to the iron-peak elements, Ca, and Ba. These re- 
sults, given in Table [6] demonstrate the overall accuracy of our 
pulsation measurements. Among the lines with definite radial ve- 
locity variation the lowest amplitude ofl5±l ms" 1 is detected in 
the strong Yll 5662.9 A line (W\ = 51 mA). In comparison, for 
the non-pulsating lines of the same intensity the amplitudes are at 
the level of 2-4 ms -1 . We consider this to be the upper limit of 
pulsation amplitude for the lower atmosphere of 10 Aql. Note that 
although for many non-pulsating lines the free-period search yields 
values comparable to real fre quencies of 10 Aql (see n ext Section), 
the False Alarm Probability dHorne & Baliu nas 1986) of such sig- 
nals is always large (Prob, listed in the Table is small). 

In addition to the centre-of-gravity RV measurements, we 
studied bisector velocities across the Ha line core and profiles of 
the strongest REE spectral lines. 

The UVES and SARG datasets overlap for about 1 h in the 
night of luly 15. The radial velocities from the two spectrographs 
agree within the error bars, where the RV measurements with the 
SARG spectra are less accurate due to a lower S/N. 



4 FREQUENCY ANALYSIS 

We carried out the frequency analysis by applying the standard 
combination of discrete Fourier transform (DFT) and least-squares 
fitting. For each of the studied REE ions we found a number of lines 
showing clear pulsation signatures. Aiming to improve the signal- 
to-noise ratio of the RV data, we subtracted linear trends from the 
measurements of different lines of the same ions and averaged these 
data. The RVs derived from the average of 5 Dy III lines turned 
out to be especially useful for the frequency analysis. In all nights 
of our UVES observations, 10 Aql exhibited a large variation in 
the RV amplitude, which often changed from < 100 m s _1 to 500- 
700 ms -1 in the course of a few hours. This behaviour, illustrated 
in Fig.[2]for Dy III, is caused by beating of several frequencies. 

The total time span of the spectroscopic data set is about 25 
hours, but the spectra were obtained in a period of 21 days and, 
thus, are characterized by a very low duty cycle. Despite the clarity 
of the pulsational RV variations, severe aliasing inhibits a secure 
identification of the frequencies from the spectroscopic data alone. 
For example, the peak with a highest amplitude is often found at 
1.3853 mHz, which is a d~ alias of the photometric frequency /j. 
Figure[T]illustrates this problem for Dy III. 

The original motivation for simultaneous MOST and spectro- 
scopic observations was to take advantage of the accurate photo- 
metric frequency information in the analysis of spectroscopic data. 
Therefore, we have adopted the three main MOST frequencies 
as the initial guess in a non-linear least-squares fitting procedure. 
Then, we optimized the amplitudes, phases and periods of the three 
frequency components for each REE ion. For a number of REE 
ions (Cell, Ndlll, Smll, EuII, Gdll, TbiII, Dylll) we were able 
to determine all three frequencies with a S/N= 10-60. The formal 
accuracy of the periods determined for each of these ions is com- 
parable or better than for the MOST photometry. We attribute this 
result to a much higher S/N achieved in the spectroscopic observa- 
tions. In Table[3]we report the weighted mean of the spectroscopic 




1.6 1.8 
Frequency [mHz] 

Figure 1. The amplitude spectra for the average Dy III radial velocities in- 
ferred from the UVES+SARG observations, a) the DFT of the radial veloc- 
ity signal; b) spectral window; c) the residual amplitude spectrum obtained 
after prewhitening the three main frequencies fi—fz (dashed lines); d) the 
residual amplitude spectrum obtained after prewhitening the four frequen- 
cies /1-/4; e) the residual noise after prewhitening the seven frequencies 
fl—ft from Table[5] These frequencies are indicated by dashed lines, while 
the dotted lines correspond to the MOST candidate frequencies. 



pulsation frequencies in 10 Aql. A comparison with MOST (Ta- 
ble [TJ shows that all our frequencies deviate by less than 1 a from 
the photometric results. 

After prewhitening the radial velocity curves with the three 
main frequencies, all REEs show residual amplitude peaks in the 
1.3-1.5 mHz range. In the UVES data the residual RV variations 
have amplitudes between 10 and 50 ms -1 , and show different fre- 
quencies at the S/N level of 3-4, depending on the ion. None of 
these frequencies coincide with the candidate MOST frequencies 
ft and fj. Accurate determination of the low-amplitude frequen- 
cies in spectroscopy is hampered by the aliasing problem. Never- 
theless, we note that the most promising spectroscopic candidate 
frequency is / = 1.3367 mHz. Variation with this period (or its 
aliases) is seen in Dy III, Ce II, Sm II, and Nd III. This frequency 
coincides with a sidelobe of a MOST orbital frequency harmonics 
and therefore does not contradict the space photometry reported by 
iHuber et"al] fe008h , because these frequency ranges were excluded 
deliberately from the analysis. The other candidate spectroscopic 
frequencies are 1.41 1-1.434 mHz (seen in Dy III, Tb III, Ce II, and 
Smll) and 1.505-1.525 mHz (Dylll, Smll, and Gdll). They are 
not observed in the MOST data. The amplitudes and phases of the 
three main periods do not change beyond their formal error bars 
when additional low-amplitude frequencies are included in the fit. 

A lower S/N and the specific temporal coverage of the SARG 
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Figure 2. The average radial velocity variation of Dy III lines during seven 
nights of spectroscopic observations. The measurements are shown with 
black (UVES) and green (SARG) dots. The solid curve represents a least- 
squares fit of seven pulsation frequencies detected in 10 Aql. 
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Figure 3. Part of the average Dy III radial velocity curve (dots) around the 
amplitude minimum. The three-frequency and seven-frequency solutions 
are shown by dashed and solid lines, respectively. 

observations close to minimum RV due to beating effects do not 
allow us to use these RV measurements for all REE ions. However, 
we found that the highest amplitude Dy III lines allow to analyse 
the complete set of UVES+SARG data and improve the frequency 
solution. 

Three practically unblended DyHI lines (5345.37, 5556.13, 
and 5730.34 A) as well as two unidentified lines (5373.02 and 
5471.42 A) that we propose to identify also as Dy III were chosen 



Table 3. Secure (bold) and tentative (italics) spectroscopic pulsation fre- 
quencies of 10 Aql determined from the analysis of the average RVs of 7 
REE ions derived from the UVES spectra (second and third columns) and 
from the average RVs of 5 Dy III lines measured in the UVES+SARG spec- 
tra (fourth and fifth columns). The numbers in brackets give the error in 
units of the last significant digit. 



id 


v (mHz) 


ion 


v (mHz) 


S/N 




UVES data 


UVES+SARG data, Dy III 


h 


1.447842(4) 


all 


1.44787(1) 


39.8 


h 


1.396883(3) 


all 


1.39684(1) 


48.2 


h 


1.427082(7) 


all 


1.42711(1) 


21.2 


h 


1.4108 


Dym,Tblll 


1.41627(3) 


7.9 




1.4341 


Cell 








1.4177 


Smll 






h 


1.3127 


Dy III 


1.35227(3) 


6.6 




1.3367 


Cell 








1.3424 


Smll 








1.3600 


Ndm 






h 






1.48260(3) 


6.0 


h 


1.5048 


Dy III 


1.50488(5) 


4.3 




1.5125 


Smll 








1.5252 


Gdll 







for RV measurements. The linear instrumental trends were sub- 
tracted for each night of observations. Because of the negligible 
phase shift between the RV curves of these lines we can aver- 
age the RVs and weight data according to the individual errors 
to increase the accuracy of the frequency analysis. After subtract- 
ing the three main frequencies, that are the same as detected by 
MOST or by using only the UVES data, the residual RV varia- 
tions show amplitudes clearly exceeding the noise level in three dis- 
tinct regions, indicating the presence of additional signal in the data 
(see Fig.[T};). To reduce the aliasing problem the program SigSpec 
teeeganll2007h was applied to the residuals for an accurate deter- 
mination of additional low-amplitude frequencies. The main func- 
tionality of SigSpec is the evaluation of the probability that a sig- 
nal in the Fourier-domain is not caused by random noise by taking 
frequency, amplitude and phase information into account, yielding 
a quantity called spectral significance. In addition, SigSpec em- 
ploys an anti-aliasing correction procedure, in which several trial 
frequencies (starting with the one showing the highest significance) 
are prewhitened. The most probable frequency is assumed to be 
that which results in the lowest residu al scatter of the data after 
prewhitening. The program Period04 dLenz&Bregei]|2004l) was 
used to check the SigSpec solution by calculating a multi-sine 
fit, allowing simultaneously to improve frequency, amplitude and 
phase of all detected signals. 

For the residual RV variations, the selection of the frequency 
at / = 1.4163 mHz yields the lowest residual scatter. Note that this 
peak does not coincide with the highest peak in the amplitude spec- 
trum (see Fig.QJ). Choosing another (e.g., the highest amplitude) 
frequency in this prewhitening step will yield an almost completely 
different set of subsequent frequency values, but not reaching the 
same low noise level. 

After prewhitening of this signal two additional power excess 
regions remain (Fig. [T]l). In total, three additional frequencies are 
needed to obtain residuals that show only noise (Fig. QJ). Table [3] 
and[5]lists all derived frequencies. 

Figure [2] shows a direct comparison of the RV observations 
and a seven-frequency solution, while Fig.[3]displays an amplified 
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Table 4. Summary of our spectroscopic pulsational analysis of individual spectral lines in 10 Aql. The columns give central wavelengths in A, followed by 
the free search period P with the respective error e stimate ap derived with a least-squares fit. The corresponding probability of a periodic signal (Prob.) 
is calculated according to iHor ne & Baliunas] jl986h . The next three groups of columns give pulsation amplitudes and phases ip with the respective errors, 
determined with a simultaneous fit of three fixed pulsation periods indicated in the column head. Phases are given as a fractional pulsation period. The starting 
point for frequency calculations is HJD=24453919. 60523. 



Free period Fixed periods 

A P= 11.931 min P= 11.511 min P= 11.679 min W\ 

(A) P (min) ap Prob. A a a f a v A a a f a v A a a <P a v (mA) 
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0.00013 
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77 


7 
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0.015 


61 


7 


0.480 


0.019 


14 


7 


0.986 


0.017 


6.3 


5044.018 


12.03110 


0.00007 


1.0000 


129 
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0.856 


0.006 


109 
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90 
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0.010 


7.1 
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12.03107 
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1.0000 


263 
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186 


5 


0.555 


0.005 


141 
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12.03110 


0.00009 


1.0000 
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0.010 


76 
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1.073 


0.013 
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5468.391 
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0.00009 


1.0000 


189 
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131 
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0.009 


100 
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1.043 


0.013 
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5512.047 


12.03094 
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1.0000 
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12.03097 


0.00016 
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56 


7 


0.650 


0.021 


36 


7 


1.197 


0.032 


23.3 


5802.534 


12.03121 


0.00010 


1.0000 


92 


4 


0.983 
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67 


4 


0.632 


0.010 


43 


4 


1.128 


0.016 


26.2 


5845.020 


12.03085 


0.00014 


1.0000 
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0.010 


52 
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14.1 
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0.00007 


1.0000 


151 
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0.007 
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6 


0.329 


0.007 


132 
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0.844 


0.009 


15.1 


5103.082 


12.03093 


0.00007 


1.0000 


145 


3 


0.692 


0.004 


146 


3 


0.323 


0.004 


107 


3 


0.836 


0.006 


16.5 


5116.686 


12.03100 


0.00008 


1.0000 


184 


8 


0.670 


0.007 


193 


8 


0.326 


0.007 


141 
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0.836 


0.010 
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5759.515 


11.42085 


0.00012 
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134 
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0.299 


0.009 


78 
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0.844 


0.017 
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1.0000 


71 
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0.027 
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0.030 


31 


12 


0.786 


0.063 


5.4 
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11.89795 
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1.0000 


63 
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0.594 


0.019 


55 


7 


0.216 


0.022 


28 


7 


0.767 


0.044 


9.3 


6303.408 


12.03093 


0.00020 


1.0000 


81 


7 


0.579 


0.015 


37 


7 


0.179 


0.034 


27 


8 


0.759 


0.046 


7.3 


6437.670 


12.03094 


0.00016 


1.0000 


72 


5 


0.588 


0.011 


48 


5 


0.226 


0.017 


32 


5 


0.818 


0.025 


26.1 


6645.070 


12.03091 


0.00009 


1.0000 


108 


3 


0.624 


0.006 


70 


3 


0.254 


0.009 


52 


3 


0.834 


0.012 


34.5 
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Table 5. Continuation of Tab.|4] 



Free period Fixed periods 

A P= 11.931 min P = 11.511 min P= 11.679 min W\ 

(A) P (min) crp Prob. A a a tp a v A a a ip a v A a a tp a v (mA) 
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11.78693 


0.00008 


1.0000 
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0.010 
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6 


0.422 


0.010 


98 


6 
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0.011 


12.2 
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1.0000 


28 


8 
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0.049 


66 
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0.021 


68 
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5.6 
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22 


1 1 
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0.082 
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1 1 


0.172 


0.018 
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1 1 
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4.7 
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1.0000 


121 


14 
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0.019 


92 


14 


0.434 


0.024 


99 


14 


0.947 
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5815.839 
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0.00022 


1.0000 


17 


6 


0.786 


0.061 


60 


6 


0.189 


0.018 


46 
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8.8 
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1.0000 


15 
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98 
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0.022 
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0.032 


5.1 
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11.42080 
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1.0000 
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0.016 
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9 
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0.005 
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0.524 


0.006 
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10 


1.030 


0.008 


9.0 


5471.416 


12.03106 


0.00009 


1.0000 


341 


13 


0.888 


0.006 
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13 


0.523 


0.008 


196 


14 


1.024 


0.012 


9.0 


5556.132 


12.03097 


0.00008 


1.0000 


349 


11 


0.875 


0.005 


311 


11 


0.537 


0.006 


208 


11 


1.045 


0.009 


7.8 


5730.335 


12.03104 


0.00009 


1.0000 


319 
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0.004 


273 
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0.005 
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1.051 


0.007 


12.0 
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5543.732 


12.03096 


0.00010 


1.0000 


77 


3 


1.098 


0.007 


63 


3 


0.750 


0.008 


39 


3 


1.278 


0.014 


11.3 



part of the observations near beating. Clearly, using only the three 
definite MOST frequencies is not sufficient to fully reconstruct the 
observed variations. While the 7 frequency fit is an improvement 
one has to stress that the frequency identification is insecure due to 
the severe aliasing effects and it strongly depends - as was already 
mentioned - on the prewhitening sequence. 

The discrepancy that none of the spectroscopically identified 
low-amplitude frequencies show up in the MOST photometry can 
be explained by their expected amplitude in the photometry of 
ab out 30ppm, which is about 3 times the noise level as stated 
bv lHuberetaIU2008t) . This estimate comes from scaling the RV 
amplitude of /1 to the highest peak in the residual RV variations 
and assuming the amplitude ratios to be similar in photometry and 
spectroscopy. Therefore, this additional signal might be lost in the 
noise. On the other hand, it has also been suggested that spectro- 
scopic and photometric observations, probing different layers of the 
stellar atmosphere, can yield different pulsation frequencies (Kurtz, 
Elkin & Mathys 2006). The clear lack of signal around 1.41 mHz 
in photometry, but which shows a convincingly high S/N in spec- 
troscopy, could support this concept. That neither of the photomet- 
ric candidate frequencies /i nor /j contribute to a good frequency 
solution for the RV variation s corroborate that t hese frequencies 
are insecure, as was stated bv lHuber et al.l (2008), or indicate that 
the amplitude scaling between photometry and spectroscopy is not 
as simple as we assumed. The spectroscopic data does not show 
any signal close to the t hird frequency 1 . 24 m Hz) proposed by 
iHeller & Kramed jl990h and lBelmonte et al.l dl99lh . This is consis- 
tent with the results obtained by MOST. 



As already noted, our dominant pulsation frequencies are con- 
sistent with those determined from the analysis of the MOST data, 
but we find systematically different amplitude ratios. The MOST 
results indicate that /1 shows the hig hest amplitude an d fa is the 
second highest: A1/A2 = 1.1 ± 0.1 l IHuber et alj|2008h . However, 
in spectroscopy fa (P — 11.931 min) exhibits the highest ampli- 
tude for nearly all lines (A1/A2 = 0.7-0.9). One can speculate 
that this discrepancy comes from a different sampling of the prop- 
agation regions in spectroscopy and photometry, respectively. 

After improving the estimate for the three high-amplitude 
pulsation frequencies in lOAql, we determined amplitudes and 
phases (expressed as fractional pulsation period) for 65 selected 
lines based on this solution and using cos {2n [(t — to ) jP — tp] } to 
fit the RV curves instead of the usual cos {27r[(t — to)/P + <fi]}, 
which correlates a larger phase to a later RV maximum. 
The reference time adopted for pulsation phase calculations is 
HJD=24453919.60523. The results of the frequency analysis of 
individual lines in the data set containing all 5 UVES nights are 
reported in Table [4] For all measured lines we also performed a 
period search with the periodogram method, which allowed us to 
estimate the probability of th e pulsation signal detecti on (1— False 
Alarm Probability (FAP), see lHorne & Baliunasll 19861) . With a few 
exceptions, the FAP is below 10 -4 . 

Tables|4]also gives the equivalent width of the pulsating lines. 
Notably, all REE lines with amplitudes significantly exceeding 
100 ms' 1 are very weak (equivalent width ~ 10 mA or less). Their 
analysis becomes possible thanks to a high S/N ratio achieved in the 
time-resolved UVES spectra of 10 Aql. 
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Table 6. Summary of the frequency analysis of non-pulsating lines in 
lOAql. The columns give central wavelengths in A, followed by the free 
search period P with the respective error estimate ap derived with a least- 
squares fit. The corresponding prob ability of the periodic signal (Prob.), is 
calculated according to iHorne & Bafiunasl ll986h . The last three columns 
give an estimate of the RV amplitude and error, and equivalent width of the 
line. 



Table 7. Phase lags (in seconds and as a fractional main period P = 11.93 
min) between the maxima of luminosity and RV variations of different 
chemical species in the atmosphere of 10 Aql. 
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5 PHASE RELATIONS BETWEEN PHOTOMETRY AND 
SPECTROSCOPY 

Our spectroscopic time-series were carried out simultaneously with 
the MOST photometry. This gives us a unique possibility to de- 
rive directly the phase lag between the photometric and spectro- 
scopic pulsational variations, which provide useful constraints for 
subsequent modelling of oscillations. Earlier attempts of similar si- 
multaneous analysis of the spectroscopic and photometric obser- 
va tions were carried out only for one roAp star - HD 24712 - 
by Ma tthews et al. ( 1988) using ground-based photometry and by 
Rvabchiko va et al.l j2007ah using MOST photometry. 

The high S/N and spectral resolution of the present observa- 
tions of 10 Aql allow us to derive phase lags for individual spectral 
lines, sampling different atmospheric layers. In order to minimize 
the influence of the higher (with respect to the spectroscopic obser- 
vations) point-to-point scatter of the photometric observations, we 
have constructed an artificial light curve based on the frequencies, 



Line 


A (A) 


seconds 


Phase lag 

fractional period 


Eu II 


6645 


40 ± 8 


U.05 ± 0.U1 


Hi 


6563 


58 ± 12 


0.06 ±0.02 


Smll 


5053 


84 ±24 


0.12 ±0.03 


Gdll 
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0.22 ±0.06 


Dy III 
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0.29 ±0.01 


Dy III 
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0.31 ±0.03 
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5274 
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0.32 ±0.02 


Dy III 
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217 ± 17 


0.32 ±0.02 


Ndm 
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260 ± 25 


0.40 ± 0.03 


Ndm 


5802 


273 ± 23 


0.40 ± 0.03 


Ndm 
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312 ± 17 


0.44 ± 0.02 


Ndm 


5294 


442 ± 30 


0.59 ±0.04 


Ndm 


5102 


444 ± 43 


0.62 ±0.06 


Prill 


5300 


449 ± 18 


0.70 ±0.03 


Tbm 


5505 


559 ± 7 


0.77 ±0.01 



amplitudes and phases of the 3 dominant pulsation modes (/i, fa 
and /j) as derived by a multi-sine fit to the MOST data. In a next 
step the artificial time-series were cross-correlated with the RV ob- 
servations of the individual spectral lines. The time interval for the 
cross-correlation was chosen from plus to minus 1 1 .93 minutes (the 
period of the main photometric frequency, /i), with an increment of 
1 second. We note that the significant modulation of the pulsation 
amplitude due to beating of several frequencies helps to distinguish 
between a phase lag of, e.g., -0.2 and +0.8, which is impossible for 
a monoperiodic star. 

The results of the cross-correlation analysis for the represen- 
tative set of spectral lines are given in Table|7] where we define the 
phase lag as the difference between the RV maximum of a given 
line and the photometric maximum, expressed in seconds or as a 
fractional main pulsation period P — 11.93 min. The error of the 
phase lag determination was estimated as follows. We added a nor- 
mally distributed random signal to both the spectroscopic and pho- 
tometric data assuming that the noise amplitude corresponds to the 
observational error. For spectroscopy the error was derived from 
the RV determinations, and for photometry the MOST team rec- 
ommended to use 2% of the signal. A phase lag was derived from 
these noisy data and repeated 200 times. The resulting standard de- 
viation was adopted as an error of the phase lag determination. 

In contrast to HD 24712, where the RV maxima for all lines 
preceded the photometric maximum, in 10 Aql the maximum of 
brightness occurs prior to the RV maxima but the phase lag itself 
depends strongly on the line considered. It is minimal (~ 40 s) for 
the Eu II 6645 A line and gradually increases for other lines, reach- 
ing a maximum of 550 s in Tb III. The growth of the phase lag from 
Eu to Tb lines has the same character as in HD 24712. It probably 
traces the outward propagation of the pulsation wave through the 
atmospheric layers. 
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6 PULSATIONS IN SPECTRAL LINES 
6.1 Phase - amplitude diagrams 

The presence of chemical stratification in the atmospheres of Ap 
stars offers a unique possibility to resolve the vertical structure 
of pulsation modes. However, the line-forming regions for REE 
are difficult to derive. Atmospheric modelling for such purposes 
should account for chemical stratification, deviations from Local 
Thermodynamical Equilibrium (LTE), magnetic field, and, eventu- 
ally, pulsation effects on the shape and intensity of spectral lines 
with substantial RV variability. The formidable difficulties of such 
an elaborated theoretical approach explain why up to now Nd and 
Pr line formation calculations accounting for NLTE and magnetic 
field effects in chemically stratified atmosphere have been carried 
out only for two roAp st ars, 7 Equ and HP 247 12 (Mas honkina, 
Ryabchikova & Rvabtsev 12003 ; IRvabchikova et alj|2007ch . 

The assumption that the line intensity or equivalent width rep- 
resents a proxy of the relative for mation heights, frequently used in 
early roAp pulsation studies (e.g. iKanaan &Hatzesll 19981) . is only 
justified for the same element, respectively ion. It is now well estab- 
lished that stratification effects are dominant in the atmospheres of 
cool Ap stars and, therefore, the formation region of the weak lines 
of one element is not necessarily located deeper in the atmosphere 
than for strong lines of another element. Consequently, discussing 
pulsation amplitude or phase versus line depth in context of pulsa- 
tion modes can be gross l y misle ading. 

iRvabchikova et all 12007b) suggested a different approach 
to the pulsation tomography problem based on the assumption 
of a continuous amplitude versus phase relation for an outward 
propagating magnetoacoustic wave. They proposed to use phase- 
amplitude diagrams to diagnose the structure of pulsation modes 
based on the following arguments: If one neglects surface chemical 
inhomogeneity (which is a reasonable assumption for slow-rotating 
roAp stars), then both phase and amplitude may provide informa- 
tion on the relative line formation depth. For instance, a phase- 
independent growth of the amplitudes measured for groups of lines 
belonging to different chemical species corresponds to a standing 
wave in the layers where these chemical species are formed. In the 
case of a phase-dependent amplitude trend, the lines showing later 
RV maximum should originate higher in the atmosphere. In the 
general case the vertical pulsation structure of roAp stars is likely 
to be a superposition of standing and running waves, correspond- 
ing to the decoupled magnetic and acous tic pulsation component s 
which have different depth dependence (Sousa &Cunhal l2008b). 
In the same star pulsations in certain ions may trace the standing, 
magnetic part of the wave while other ions may show predomi- 
nantly the running, acoustic component. The phase-amplitude dia- 
gram method is well-suited to characterize this behaviour thus of- 
fering the possibility of tracing the vertical properties of pulsation 
waves and studying the physics of magneto-acoustic oscillations 
from trends of pulsation velocity amplitude and pulsation phase 
without tedious assignment of the physical depth to each pulsation 
measurement. 

The key assumption used in interpreting the phase-amplitude 
diagrams is that a later RV maximum corresponds to a higher 
line formation. This is supported by the gradual increase of pul- 
sation phase with height from H/3-, H q-cores to Ndn-III - 
Pr II- III lines (e.g.. ISachkov et alj|2006h in practically all roAp 
stars where formation heights of pulsating lin es could be con- 
strain e d independently of the puls ation analysis ( Mash onkina et al.l 
120051 : IRvabchikova et alj |2007c). Since the chemical proper- 
ties of most roAp stars are quite similar tevabchikova et all 



2004) and their phase-amplitu de dependences are also compara- 
ble I Ryabchikova et al]|2007bl) . the assumption of continuous out- 
ward increase of pulsation phase appears to be reasonable from 
the empirical point of view. On the other hand, the general va- 
lid ity of this picture has b een questioned in the theoretical work 
bv lSousa & Cunhal |2008a). These authors demonstrated that using 
a certain combination of the model parameters (not applicable to 
lOAql or to any other roAp star), one can obtain a complex su- 
perposition of standing and outwardly propagating waves which 
would mimic an inwardly running wave. This situation may com- 
plicate interpretation of the phase-amplitude diagrams, and the case 
of 33 Lib (Sect. |9.3T l possibly represents an example of such com- 
plex situation. However, a more realistic theoretical modelling us- 
ing pulsation and atmospheric parameters relevant for specific roAp 
stars is required to confirm the reality of such complex depth de- 
pendence of pulsation phase. 

Similar to the analy sis of 10 roAp stars presented by 
IRvabchikova et all (12007b). we applied the amplitude-phase dia- 
gram method to lOAql. Furthermore, the availability of our rich 
spectroscopic observational material allows us to improve the pre- 
cision of our previou s analysis which was presented at the CP#Ap 
Workshop in 2007 bv lSachkov et alj ( 120071) . and for the first time to 
recover the amplitude-phase diagrams for several pulsation modes 
in the same roAp star. These diagrams are shown in Fig.|4]for the 3 
main frequencies observed in 10 Aql. 

For the purpose of comparison, we plot bisector measurements 
across the Ha core in the amplitude-phase diagrams constructed 
via centre-of-gravity measurements. These measurements provide 
us with a physical depth scale in the atmosphere because hydrogen 
is distributed homogeneously with height. The core of the Ha line 
is formed in the region between logrsooo ~ —2 and logrsooo ~ 
—4.7 according to NLT E calculations using the codes described in 
iMashonkina et al.l J2007l) . 

As in other roAp stars, measurable pulsation amplitudes ap- 
pear in the layers where Eu II, La II, and the core of Ha are formed, 
reach maximum and then show decrease of the amplitude. The 
phase increases by 0.2 across the Ha-line core while the RV am- 
plitude grows by 4-5 times. We observe a running wave-like be- 
haviour in this part of the stellar atmosphere. Above the formation 
height of the deepest point of the Ha profile (r = 0.15) the am- 
plitude continues to grow rapidly without significant phase change. 
However, while in most other roAp stars the maximum RV ampli- 
tudes are observed in singly and doubly ionised Nd and Pr lines, 
in 10 Aql the maximum are seen for the lines of Ce II, Dy III and a 
couple of unidentified lines too, attributed by us to Dy III. 

The lines of Nd and Pr are formed above the Dy III lines 
and show a RV amplitude rapidly falling to very small values. In 
even higher layers, the pulsation amplitude starts to grow again 
and reaches a new maximum in lines of Tb III. Simultaneously, 
the pulsation phase changes by ~0.5 between the two maxima. 
This is a typical picture of a pulsation node in the stellar atmo- 
sphere: the phase changes by 180° while on both sides of the 
phase jump we observe comparable RV amplitudes. The presence 
of the node is supported by the bisector measurements in Ndlll 
lines (see below). If the optical depth scale indicated by the core 
of the Ha line is valid for other lines too, then we expect that 
the pulsation node is located near or above the optical depth of 
logrsooo ~ —4.5. However, it should be stressed that the ampli- 
tude and phase distributions with height in the atmospheres of roAp 
stars is exp ected to depend on the lo cal magnetic field strength and 
inclination l lSousa &Cunhal l2008b). Thus, the inferred location of 
the node observed in the disk-integrated quantities represents the 
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4. Amplitude-phase diagrams for three main pulsation periods in 
P = 11.93 min (a), P = 11.51 min (b), and P = 11.68 min (c). 



average of the local vertical pulsation structure which, possibly, 
exhibits substantial variation over the stellar surface. Ideally one 
should take this variation into account when determining the exact 
vertical location of the node. 

Among the lines of elements other than rare earths only the 
strongest line of Y II, A 5662.9 A, show s low-amplitude but definit e 
pulsation. As in most other roAp stars tevabchikova et al]l2007bl) . 
this line overlaps with Pr III lines in the amplitude-phase diagrams. 

The amplitude-phase diagrams for the three pulsation modes 
are very similar except, perhaps, Gdll and La II lines. Only 4 out 
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Figure 5. Bisector RV amplitudes (top) and phases (bottom) across the 
Ha core for the three main pulsation periods in 10 Aql: Pi = 11.51 min 
(crosses), P2 = 11.93 min (filled circles), and P3 = 11.68 min (dia- 
monds). The upper x-axis indicates the optical depth scale. 



of 14 Gd II lines have equivalent widths larger than 7 mA. Also, 
most of Gd II lines are shallow due to substantial magnetic splitting; 
their central depth does not exceed 4-5 %. This makes pulsation 
measurements difficult. Pulsations can be detected in all Gd II, and 
it seems that this element shows higher amplitudes for the 11.51 
and 1 1.68 min periods. The same situation is found for La II lines, 
which have even lower equivalent widths and lower RV amplitudes. 
This discrepant variation of the Gd II and La II lines for different 
frequencies may indicate a difference in the vertical extension of 
the corresponding modes. In particular, on the basis of these mea- 
surements, one may suspect that the frequency fa (P = 11.93 min) 
corresponds to a pulsation mode located in higher atmospheric lay- 
ers compared to fa and fa. This agrees with the difference between 
photometric and spectroscopic amplitude ratios for fa and fa. 

6.2 Bisector variation 

The RV amplitudes and phases across the Ha core are shown in 
Fig.|5]as a function of the normalized flux (scale at the bottom) and 
optical depth (scale at the top) for all three principal pulsation pe- 
riods seen in 10 Aql. We find very similar behaviour for the three 
modes. In fact, the amplitudes of fa and fa are identical within 
error bars, whereas the amplitude of fa is slightly lower. This sug- 
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Figure 6. The RV amplitudes (top) and phases (bottom, expressed in as a 
fractional pulsation period) for Nd III 5 102 A. 

gests that the three modes sample the hydrogen line core formation 
region in a similar manner. 

Because of the weakness of most spectral lines, the only metal 
line suitable for a precise bisector analysis is Nd III 5 102.43 A. This 
line is deep enough and its wings are free from blends. In many 
roAp stars this Ndm line is blended with Ndll 5102.39 A and 
should be used with caution. However, due to the overall weakness 
of the Nd II spectrum in 10 Aql, the blending effect is negligible. 

Bisector measurements across the Ndm 5102 A line profile 
are shown in Fig. [6] A phase jump of approximately 0.4-0.5 seems 
to be located in similar atmospheric layers for all three main fre- 
quencies. A bisector analysis provides an independent evidence for 
the pulsation node in the u pper atmosphere of 10 Aql. This is the 
second case after 33 Lib jMkrtichian et alj 2003; Kurtz, Elkin & 
Mathy s 120051 : iRvabchikova et al J 2 007b ) of a pulsation node found 
in a roAp star located within the REE-rich cloud. 

Ndm 5102 A line provides the most clear signatures of the 
pulsation node: phase jump of ~0.5 of the period accompanied by 
the minimum in RV distribution. The start of the phase jump is s een 
in other Nd III lines as illustrated by Fig. 6 of Elk in et alj J2008h . In 
Nd III 6145 A line it also reaches ~ 180° , however, we did not an- 
alyze this line due to its known blend with Sii line that influences 
bisector RV values. Another strong Ndm 6237 A line was omit- 
ted from the bisector analysis, too, because it is blended in the red 
wing by one of Sm II, an element which also shows pulsation, but 
probably is formed deeper in the atmosphere (see our Fig.0. 

In comparison to the bisector variation in the Ha core, the 
Ndm 5102 A bisectors behave differently for the three modes. 
The amplitude in the outer wings is the highest for /2, whereas 
the steepest phase jump is seen for fa. Based on the location of the 
maximum change of the bisector phase further outside in the wings, 
we tentatively conclude that the node region is located somewhat 
deeper for the fa mode than for the other two frequencies. 



7 ND III LINE PROFILE ASYMMETRY 

Many cool Ap stars exhibit an unusually large broadening of the 
doubly ionized Nd and Pr lines, for which macroturbulent veloc- 
ities of the order of 10 kms -1 are often needed to achieve a 
good fit to the observed line profiles iKochukhov & Rvabchikovi] 



l2001al ;l Rvabchikova et al ] |2007bl) . The nature of this phenomenon 
is poorly understood, but it could be related to a turbulent instabil- 
ity in the upper atmospheric layers, possibly caused by the anoma- 
lous temperature gradient produced by a stratified distribution of 
the chemical elements. Although roAp stars are the ones most fre- 
quently showing this anomaly, it does not correlate with pulsa- 
tion amplitude and is equally strong in stars in which pulsations 
at the level of 20-30 m s -1 are barely detectable (Kochukho v et alj 
2008). As discussed bv lKochukhov et alj J2007I) the modulation of 
turbulence by the periodic contraction and expansion of the pulsat- 
ing atmosphere in roAp stars with higher amplitudes explains the 
characteristic a symmetric line profile variat i on of t he REE lines, 
first detected by Kochukho v & Rvabchikovj d2001ah in 7 Equ. 

10 Aql, too, shows remarkably broad lines of Ndm, but 
curiously we find that the strongest of these lines are also 
extremely asymmetric. This effect is illustrated in Fig. [7] for 
several Ndm lines. In this figure we show rare-earth line 
profiles in the spectra of four slowly rotating roAp stars: 
two stars with a small magnetic field, 10 Aql (X' e ff = 7550 K, 
logg = 4.0, 7j c sini = 2.0 kms -1 , (S) = 1.5 kG) and HD 12932 
(T eff = 7620 K, log# = 4.2, « c sini = 2.0 kms" 1 , (B) = 1.7 kG), 
and two stars with a large magnetic field, 33 Lib (T e g = 7550 K, 
log 5 = 4.3, i> c sini^2.0 kms -1 , {B}=5 kG) and 7 Equ 
(r eff = 7700 K, log q = 4.2 w e sini<1.0 kms" 1 , (B)=4A kG) 
( Rvabchiko va et all [2007b). For a comparison the profile of 
the magnetically insensitive non-pulsating (or weakly-pulsating) 
Fei 5434 line is also shown for all stars. The wavelength scale is ex- 
pressed in velocity units relative to the line centres. The strongest 
Ndm lines in 10 Aql at A 5102 and 5294 A have a remarkable 
blue asymmetric wing reaching out to —50 kms" 1 . The red wing 
stops at ~20 kms -1 as do both wings in the Ndm lines in the 
other roAp stars shown in Fig.[Vj In 10 Aql the blue wing asymme- 
try is still prominent in Nd III 6145 A and gradually disappears for 
weaker Nd III lines. The profiles of strongly and weakly pulsating 
lines in HD 12932 and 7 Equ are symmetric. In 33 Lib the weakly 
pulsating Fei 5434 A line is symmetric while REE lines have a 
small red asymmetry. In 10 Aql the non-pulsating Fei 5434 A line 
has a clearly asymmetric red wing. 

The available information is insufficient to attribute REE line 
asymmetries either to non-radial pulsations or to an anisotropic tur- 
bulent velocity field. The sign of the asymmetry suggests an out- 
ward flow of Nd-rich material in the layers probed by the wings 
of strong Ndm lines. However, this effect can be mimicked by a 
nonisotropic velocity distribution of an i on stably stratified in the 
upper stellar atmosphere, as discussed bv lMichaudl ( 119781 ) for Sr in 
hotter CP stars. A similar investigation is needed for REE in cool 
Ap stars to shed light on the abnormal shapes of strong Nd lines in 
10 Aql. 



8 ROTATIONAL MODULATION 

A comparison of the average spectra obtained in five UVES nights 
over a period of about 1 month did not show any hint of line profile 
or radial velocity changes on time scales of days. This supports the 
study by Rvabchikova et alj J2005h who have used own mean lon- 
gitudinal field measurements and data taken from the literature to 
argue for a long rotational period of 10 Aql. To search for long-term 
variations we have compared the average UVES spectra with obser- 
vations obtained in 2001 with the Gecko coude spectrograph at th e 
3.6 m Canada-France-Hawaii telescope ( Kochu khov et al J l2002h . 
We introduced additional broadening to the latter spectra to com- 
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Figure 7. REE and Fei line profiles in the spectra of the weakly magnetic roAp stars 10 Aql, and HD 12923, and of stronger magnetic roAp stars 33 Lib, and 
7 Equ. 



pensate for the difference in spectral resolution between Gecko 
and UVES. The global average spectrum of 10 Aql in the 6140- 
6165 A region and the residual spectra to this average are presented 
in Fig.Hl 

A signature of weak variability can be seen in the centres of 
few strong lines of Cat, Fe II, Ba II, Nd III, as well as in the moder- 
ately weak Cr II line. We are confident that this variability is not an 
artifact due to using different spectrographs, because the intensity 
variations of Cat, Fe II, Ba II lines are in antiphase with the varia- 
tions of Nd III and Cr II lines. Possibly, this is another evidence for 
a rotational modulation in the spectra of 10 Aql, indicating a long 
rotation period. 



9 DISCUSSION 

9.1 Reassessment of asteroseismic models 

The additional spectroscopic frequencies could offer a new view 
on the large frequency separation Av, a crucial factor for astero- 
seismology which is directly connected to the mean density in the 
star and describes the separation of consecutive radial overtones 
for high-order acoustic pulsation. Figure [9] shows a schematic am- 
plitude spectrum including all intrinsic and candidate frequencies 
of both MOST photometry and spectroscopy. To first view, there is 
no apparent equal spacing visible. If we consider the four frequen- 
cies with high S/N (solid lines) to estimate the spacing, the only 
reasonable solut ion corresponds to /i — fe = 51 /iHz, a value al- 
ready noted by [Huber etal] J2008h . Lower values (such as 20 or 
30 /iHz) although seem tempting, would however contradict the 
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Figure 9. The observed secure (black lines) and tentative frequencies from 
the photometry (blue lines) and spectroscopy (red lines). The expected posi- 
tion of radial orders for three different spherical degrees I (vertical dashed, 
dotted and dashed-dotted lines, respectively) is calculated based on a large 
separation of /i — f% = 51 /iHz. 



previou sly determined tem perature and luminosity of lOAql (see 
Fig. 4 in lHuber et al.l2008h . Assuming Au = 51 /iHz to be correct, 
Fig.|9] shows the expected position of other radial orders for three 
different spherical degrees I (vertical dashed, dotted and dashed- 
dotted lines, respectively) based on the position of the four highest 
S/N peaks. Evidently, the agreement is not very convincing. Nev- 
ertheless, an average deviation of about 5 /iHz (which is needed 
to align the observed values to the expected position) is about the 
order of the frequency shifts of consecutive radial overtones due 
to the magnetic field perturbation as predicted by theory (see, e.g., 
ICunhj|2006l) . It must be noted, however, that such suggestion can 
lead to (almost) any desired result and therefore have to be consid- 
ered with extreme caution. 

Generally, the frequencies we derived from spectroscopy do 
not contradict to those from MOST photometry. But spectroscopic 
data provi de an additional i nformation useful for modelling. Ac- 
cor ding to lHuber et al] (2008) the best model describing the pulsa- 
tion features in 10 Aql derived from spectroscopy corresponds to a 
star of 1.95M0, T e g = 7730 K and normal solar composition with 
the boundary reflection layer at log T5000 ~ —4. The spectroscopic 
analysis reveals a pulsation node in the Nd II - Nd III line formation 
region. A Non-LTE analysis of these lines suggested that they are 
formed in roAp stars in atmospheric layers above logrsooo = — 4 
( Mashonkin aet al.ll2005h . thus indicating the position of the node. 
Indeed, a pulsation node is predicted at log T5000 == —4 in a star of 
1.95M0, T c ff =7730 K but assuming helium depletion (Saio, pri- 
vate communication), which would better correspond to the typical 
chemistry of a magnetic peculiar star. Therefore, combining simul- 
taneous photometric and spectroscopic observations, we expect to 
improve the pulsation model of 10 Aql. This analysis will be pre- 
sented in the next paper. 

9.2 Interpretation of spectroscopic pulsation measurements 

Here, as in many other recent studies of pulsational RV variations of 
roAp stars, we interpret the measurements of amplitude and phase 
in terms of outward propagation of pulsation waves in a chemi- 
cally stratified stellar atmosphere. Previously, this interpretation of 



spectroscopic observations of multiperiodic roAp stars was hand- 
icapped by short spectroscopic time series which did not allow 
to resolve frequency patterns of multiperiodic pulsators. Here we 
overcome this difficulty by combining time-resolved observations 
acquired during several nights, but with a poor duty cycle, with 
a continuous photometric monitoring by MOST at the same time. 
This strategy allows us to resolve the vertical structure of the three 
principal modes in 10 Aql and to study the propagation of pulsation 
waves independently for all three frequencies. 

We detect a radial node in the region sampled by the formation 
heights of Nd and Pr for all three frequencies. The maximum am- 
plitudes are observed for the Dy III and Tb III lines that are presum- 
ably formed below and above the region of the node. The existence 
of a node is supported by the bisector analysis of the Nd III 5 102 A 
line. 

By attributing the differences in the RV curves of different 
REE ions to the vertical structure of oscillations we have implicitly 
assumed that different species sample the horizontal structure of 
pulsation modes in a similar manner. This might not be the case 
if the vertical distribution of pulsation amplitude and phase de- 
pends strongly on the stellar surface position due to a spotty el- 
ement distribution combined with an intrinsic dependency of the 
mode structures on the orientation and strength of the magnetic 
field. The importance of the latter effect, arising from a superposi- 
tion of magnetic and acousti c components of pulsatio n waves, has 
been recently emphasized bv lSousa & Cunhj J2008bl) . 

Due to a long rotation period we are unable to constrain the 
surface abundance distribution of 1 Aq l in the same way a s it could 
be done for rapidly rotating roAp stars jKochukhovll2006l) . For this 
reason it is not straightforward to distinguish between the verti- 
cal and horizontal structural effects. At the same time we see only 
small changes in the line profiles in a comparison of observations 
obtained in 2001 and 5 years later. Furthermore, the independent 
evidence for a node provided by phase-amplitude diagrams and bi- 
sector analyses of different REE ions suggests that the observed 
phase variation corresponds rather to a vertical structure and not to 
a different horizontal sampling of the stellar surface. 

9.3 Comparison with 33 Lib 

A comparison of the pulsation properties of 33 Lib and 10 Aql 
is of considerable interest because these are the only roAp stars 
with clear signatur e s of a radial node in the upper atmosphere. 
Rvabchikov a et al] i2007bl ) performed a pulsation analysis of 
33 Lib using phase-amplitude diagrams, which enable a direct com- 
parison of the vertical pulsation structure in two stars. The top panel 
of FigJjJJ displays the phase-amplitude diagram for 33 Lib, while 
bisector measurements of Nd III lines are compared for both stars 
in the bottom panel of Fig.UOl 

A pulsation node is clearly present in the atmospheres of both 
stars but the phase jump has opposite sign. Taking into account that 
the REE elements are concentrated in the upper atmospheric layers 
close to or even above the Ha line core formation zone, we can 
conclude from Figs.|4]and[T0]that the phase changes by 0.5 between 
the formation zones of Dy III and Tb III in the upper atmosphere of 
10 Aql, while it jumps by —0.5 between Nd II and Nd III in 33 Lib. 

Interestingly, the bisector measurements give us just the op- 
posite phase variation with depth if one treats the line profile as 
being produced in a normal stellar atmosphere where the line core 
is formed higher than the line wings. As we show in the lower panel 
of Fig. [TO] the phase decreases from the wing to the core in 10 Aql, 
and increases in 33 Lib. At present we do not have a definite in- 
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terpretation of this discrepant behaviour. One can suspect that the 
the broad wings of strong REE lines actually originate in the outer 
atmosphere and hence show pulsational characteristics of these lay- 
ers. 

Although the atmospheric parameters of 10 Aql and 33 Lib are 
similar, 33 Lib has a much stronger magnetic field, (B) = 5.0 kG, 
and it has a higher overabundance of REEs (see Rvabchiko va et al.l 
|2004|) . In addition, it has a shorter main pulsation period, P = 
8.27 min, with the first harmonic exhibiting the highest amplitude 
close to the position of a radial node. The negative phase jump and 
the shape of the amplitude-phase diagram (Fig. [10] upper panel) 
may be interpreted as a superposition of standing and running pul- 
sa tion waves mimicing an i nwardly propagating wave, as discussed 
bv lSousa & Cunh a ( 2008b). A detailed study of chemical stratifica- 
tion and atmospheric structure of both stars is required for a secure 
interpretation of pulsation results and subsequent theoretical mod- 
elling. 



9.4 Comparison with previous spectroscopic pulsation 
studies of 10 Aql 

Non-radial pulsations in 10 Aql were investigated in four differ- 
ent spectroscopic observing campaigns, including the present one. 



iKochukhov et al.1 d2002h found that RV amplitudes in EuII and 
Gd II lines exceed those determined f or Nd III line, which was not 
typical for a roAp star. Additionally, IKochukhov et al.l d2002h de- 
tected a change of RV amplitude between two consecutive nights 
of their observations. These characteristics of pulsational behaviour 
of 1 Aql are confirmed in the pr esent more extensive analysis. 

Hatze s~& Mkrtichianl |2005) observed the star during 3 nights, 
clearly detecting pulsational variability in 5 lines. They found the 
highest pulsation amplitudes of a few hundred ms _1 in the three 
spectral lines of AA 5373, 5471, 5730 A, which could not be iden- 
tified. We also find significant amplitudes in these spectral features 
and tentatively identify these lines as Dy III, which allowed us to in- 
vestigate variation of these lines in context of pulsation wave prop- 
agation through a stell ar atmosphe r e with chemical stratification. 

A recent study bv lElkinet ail [2008) presented measurements 
of the pulsation amplitudes and phases based on UVES spectro- 
scopic time-series over two hours, which we included in our present 
analysis allowing for an increase of our spectroscopic data by 10%. 
lElkin et aT] d2008h assumed that 10 Aql pulsates with a main fre- 
quency of 1.428 mHz, which is, in fact, close to one of the three 
main frequencies unambiguously detected in our spectroscopy and 
in the MOST photometry, but with the lowe st amplitude. The sec- 
ond frequency detected bv lElkin et alj |2008J) at 1 . 309 mHz is based 
on variations of Tb III lines, and is not seen in our more voluminous 
data. It probably is an artifact caused by, e.g., the limited time span 
of their spectroscopic observations of this multiperiodic pulsator. 

Our study is the first to resolve the frequency spectrum of 
10 Aql with RV measurements. We find, in agreement with precise 
space photometry, that pulsations of 10 Aql are dominated by three 
modes of comparable amplitude. Based on simulteneous spectro- 
scopic and photometric observations, we were able to infer pulsa- 
tional characteristics of individual modes as a function of height in 
the atmosphere of 10 Aql. We also obtained phase shifts between 
luminosity and RV variations providing constraints for modelling 
of roAp pulsation. 

Since our time-series analysis is based on a set of three fre- 
quencies, the RV amplitudes that we report for individual lines 
cannot be directly compared with the results of previous studies 
which assumed monoperiodic v ariation of the star. However, a good 
agreement with the results by lElkin et alj ( 120081) is found if we 
compare relative amplitudes and phase differences of the spectral 
features in common. The two studies obtained qualitatively simi- 
lar phase-amplitude diagrams, with the excepti on of discordant be - 
haviour of weak ErIII 5903 A line studied by lElkin et all d2008l) . 
The pulsation characteristics of this line deviate from the generally 
smooth phase-amplitude behaviour of other REE ions. We omit- 
ted this line from our pulsation analysis because it is located in the 
spectral region full of telluric lines and is blended by 3 atmospheric 
lines of different intensity. This blending is the most likely the rea- 
son for the deviating behaviour of the spectral line Er III 5903 A. 

9.5 Amplitude modulation in spectroscopy 

Time-resolved spectroscopic observations of roAp stars seemed to 
indicate that the pulsation amplitudes of REE lines can be mod- 
ulated on a time scale of a few h ours. This was first noted by 
IKochukhov & Rvabchikov3 d2001b|) for aCir and later observed 
by iKurtz et al.1 ((2006) for a few other roAp stars. Since in many 
cases this modulat ion cannot be link ed to known photometric pul- 
sation frequencies. lKurtz et a l. (2006) speculated that the amplitude 
modulation in spectroscopy suggests the discovery of a new type 
of pulsational behaviour in the upper atmospheres of roAp stars. 
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However, using relatively short and typically 2-hour long time se- 
ries resulting in only 100-150 individual measurements, this dis- 
covery is doubtful, because such short observational data sets do 
not allow to resolve the frequency spectrum of multiperiodic roAp 
stars. In this context it co mes not as a surprise that for the only 
two roAp stars, HD24712 1 Rv abchikova et alj|20073) and lOAql 
(this paper), for which an extensive spectroscopic (and photomet- 
ric) monitoring over many nights was performed, no unexplained 
spectroscopic amplitude modulation could be found. All amplitude 
modulations in 10 Aql can be explained by beating effects of close 
frequencies. Hence, we recommend extreme caution in the inter- 
pretation of amplitude changes seen in short data sets. 
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